Gate-Controlled Spin-Orbit Quantum Interference Effects in Lateral Transport 
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In situ control of spin-orbit coupling in coherent transport using a clean GaAs/AlGaAs two- 
dimensional electron gas is realized, leading to a gate-tunable crossover from weak localization to 
antilocalization. The necessary theory of 2D magnetotransport in the presence of spin-orbit coupling 
beyond the diffusive approximation is developed and used to analyze experimental data. With this 
theory the Rashba contribution and linear and cubic Dresselhaus contributions to spin-orbit coupling 
are separately estimated, allowing the angular dependence of spin-orbit precession to be extracted 
at various gate voltages. 



An important component along the path toward realiz- 
ing quantum "spintronic" devices 0,0] is a structure that 
allows manipulation of electron spin without destroying 
phase coherence. Spin-orbit (SO) coupling has been the 
focus of recent studies because of its potentially useful 
role in coherent spin rotators 0] , spin interference devices 
, and spin-filters 0] . The mechanisms by which SO 
coupling affects transport E3 have recently been 
considered in the context of Aharonov-Bohm (AB) phase 
and Berry phase 0, EH El El 13 EE El , underscoring 
the richness of the underlying physics. 

The conductivity of low-dimensional systems shows 
signatures of quantum interference that depend on mag- 
netic field and SO coupling El El El In par- 
ticular, constructive (destructive) backscattering associ- 
ated with pairs of time-reversed closed-loop electron tra- 
jectories in the absence (presence) of significant SO in- 
teraction leads to negative (positive) magnetoresistance 
effects known as weak localization (antilocalization) 0]. 
In this Letter, we demonstrate in situ control of SO cou- 
pling in a moderately high mobility GaAs/AlGaAs two- 
dimensional electron gas (2DEG), inducing a crossover 
from weak localization (WL) to antilocalization (AL) as 
a function of an applied top-gate voltage (see Fig. 1). 
Theory beyond the diffusive approximation must be used 
to extract gate-voltage-dependent SO parameters from 
magnetotransport when the SO precession frequency be- 
comes comparable to the inverse transport scattering 
time (t _1 ) as occurs here, and when the magnetic length 
becomes comparable to the mean free path. Theory that 
accounts for AB-like spin phases and spin-relaxation [2(j 
is developed here and used to estimate separately the var- 
ious SO terms (Rashba, linear and cubic Dresselhaus) 
over a range of gate voltages, ranging from WL to AL. 
Conventional WL theories assume SO times much longer 
than r H and so cannot be applied to clean ma- 
terials. Previous theories that go beyond the diffusive 
approximation do not treat SO [2ll l22|. or treat it only 
as spin-relaxation jH 0] without accounting for Berry 




FIG. 1: (a) Experimental magnetoconductance, Act = o(B)— ct(0), 
(circles) offset for clarity, along with three-parameter fits to Eq. 
(solid curves) for several gate voltages. Inset: Experimental mag- 
netoconductance data for the most negative gate voltage, showing 
pure WL. (b) Density and mobility as a function of V g , extracted 
from longitudinal and Hall voltage measurements, (c) Experimen- 
tal conductivity, showing strong dependence on V g . Note that 
Act ~ 1CT 3 ct. 



phase effects which play a crucial role, as we show here. 
Previous experiments in which SO rates are measured 
using WL/AL in a gated GaAs heterostructure have not 
reported in situ gate control 0, |H |2(| . Koga et at 
[27| demonstrated gate controlled SO coupling in InGaAs 
heterostructures using WL/AL. Modification of Rashba 
SO coupling using gated quantum wells has been ob- 
served using beating pat terns in Shubnikov-de Haas os- 
cillations in InGaAs [Hill, InAs/AlSb HI and HgTe 
[3l| . Gate controlled SO coupling in GaAs 2D hole sys- 
tems [H IH Q has also been investigated using beating 



2 



of Shubnikov-de Haas oscillations. The angular variation 
of SO coupling in GaAs quantum wells has been mea- 
sured using Raman scattering [35j| . 

The Hamiltonian for conduction band electrons in a 
[001] 2DEG is H = + (a ■ il), where m* is the ef- 
fective mass, k = |k| (k = (k x ,k y )) is the in-plane wave 
vector, a — (cr x ,a y ,a z ) is the Pauli spin operator and 
f2 = (Cl X) Cly) is the total SO frequency, l~i can be writ- 
ten as the vector sum of linear (Jim) and cubic (fi.D3) 
Dresselhaus terms and the Rashba term (J1_r), 



&di = oci (—xk x + yk y ) /K, 
fi R = a 2 (xfc y - yk x )/h, 
fl D3 = j(Stk a 



fc y - yk x k y )/ h - 



(la) 
(lb) 
(1c) 



where 7 arises from the lack of inversion symmetry 
of the GaAs crystal, while cti = 7(&f) also depends 
on the thickness of the wave function in the quantiza- 
tion direction, a 2 depends on the potential profile of 
the heterointerface. We assume the effect of gate volt- 
age, V g , on n (= I fl |) is through the carrier density, 
n = k 2 /2ir. Previous studies of SO coupling in single- 
interface heterostructures |3(j| support this assumption. 
The magnitude of a 2 in a single-interface heterostruc- 
ture originates mainly from the band-offset at the het- 
erointerface, which is roughly independent of V g [3fll3^|. 
The symmetry of the linear (in k) SO terms, SIdi and 
SIr, allows these terms to be represented as a spin- 
dependent vector potential A that affects the orbital mo- 
tion a nd p hase of electrons, cr • (SIdi + Or) oc k • A 
H El El El El E3 That is, the linear terms affect 
electronic interference as a spin-dependent AB-like ef- 
fect. In contrast, the cubic term, Eq. (lc), upon remov- 
ing terms with the symmetry of Eq. (la), only causes 
spin relaxation in the diffusive regime (although it also 
can produce AB-like effects in the quasi-ballistic regime 
0). To develop the theory of 2D magnetotransport 
with SO coupling beyond the diffusive approximation 
H3, we follow Refs. [HEIE ^ which treat the quasi- 
ballistic case £b < £ (Ib = \/h/2eB is the magnetic 
length and i is the transport mean free path) without 
SO coupling. The approach is to introduce an operator 
P = Gf +U (ri, r 2 , (Ti)G^(ri, r 2 , a 2 )h/2irvT for the proba- 
bility of an electron to propagate both forward and back- 
ward along a path segment from ri to r 2 , where G R (G A ) 
are single-electron retarded (advanced) Green functions, 
u 1(2) are the Pauli spin operators for particle moving 
forward (backward), v is the density of states per spin, 
and t is the scattering time. The interference contribu- 
tion from traversing a closed trajectory with n scattering 
events is given by the trace of (P) n ■ In the presence of 
SO coupling, Eq. (JIJ, the formulas in remain valid 
once a summation over spins is included in the trace. 
Introducing the total spin of interfering particle waves, 
S = (Ji + cr 2 , we write Tr[{P) n ] = \Tr[(Px) n - (P ) n \, 
where operators Pq and Pi describe singlet (S = 0) and 



triplet (S = 1) contributions. To calculate Tr[(Pom) n ], 
we diagonalize Pq(i)- We find that when Jim and Qr 
are taken into account, Pqci) has the same eigenfunctions 
as the Hamiltonian TL for particles with charge 2e, spin S 
and spin frequency 2fi: TL = ^fO 5 - - 2eA em + 2eA s ) 2 , 
where A em is the vector potential associated with the 
applied perpendicular magnetic field, B, and A5 = 
5jrp-(— &iS x — otzSy, a.2S x + otiSy) is the SO vector po- 
tential. For 5 = 0, the eigenstates are Landau levels for 
a charge 2e particle in the magnetic field B, analogous 
to the spinless problem J23|. For 5=1, eigenstates of TL 
and Pi in general require a numerical solution, although 
analytic solutions exist when either ct\ or a 2 equals zero 
[20| . An analytic solution is also found for an, a 2 7^ 0, 
when Ib < A so , where X so = (2a 1 ( 2 )m*/fi 2 ) -1 is the dis- 
tance over which spin rotates appreciably (if £ > A so ) or 
dephases (if £ < \ so ) due to spin AB-like effects. Per- 
forming a unitary transformation TL — > TL = U^TLU , 
with U = exp (— ieAs ■ r), and expanding in coordi- 
nates, we find TL = j^r(k — 2eA em + S z &) 2 , where 
a = H cS rxz/(2h 2 ), and H cS = 2(a\- a\)m* 2 j eh z is the 
effective SO field. P\ can then be block-diagonalized for 
each m (m = 0, ±1) using the Landau basis for particles 
with charge 2e in the magnetic field B — mH c g. Thus, the 
effect of fini and Hp is to produce AB-like spin phases 
El El El El E3 Higher expansion terms to TL de- 
scribe spin flip processes and can be taken into account 
by introducing a spin relaxation time r so and its corre- 
sponding field scale H so = fvr/ (2e£ 2 T so ). The resulting 
quantum interference contribution takes the form 20] 



Aa(B) 



Air 2 h 



22 C(x lm ,f lm )-C(x 0Q ,f 0Q ) 



(2) 

where xs m = (B — mH cff )/H tr describes the AB dephas- 



ing in H eS , C(x,f Sr , 
V Jo°° exp(-yf Sm t - 



PN(fsm) 



,N=0 l-P N (f Sm ) ' 

t 2 /2)LN(t 2 )dt, Ln(z) are Laguerre 
polynomials, y = (2/lxl) 1 / 2 , and H tr = h/(2e£ 2 ). The 
dephasing factors fs m are given by /i±i = (1 + (H v + 
H S o)/H tr ); foo = (1 + H v /H tr )\ fio = (1 + [H v + 
2H so) I Htr), where H 9 = ft/(4eL 2 ) and L v is the phase 
breaking length. Equation J2J does not include all in- 
dependent interference terms, notably excluding Cooper- 
channel contributions due to electron-electron interac- 
tions |l8| and a reduction of WL due to electron diffrac- 
tion effects 21]. Also, in an attempt to capture the ef- 
fects of cubic terms on H e g and H so , we introduce an 
effective vector potential A * s — A5 + 7^2 (k y 2 , —k x 2 ) ~ 
A5 + 7^-j-(fc 2 , — k 2 ) which leads to an effective SO field, 

H* B = 2{a\ -a\ + 2nn ail - n 2 j 2 n 2 )m* 2 /eh 3 . (3) 

Equation © is applicable when B > H* s (see Fig 2). We 
have confirmed that fitting only to data where B > H* s 
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FIG. 2: Spin-orbit effective fields, H* a (filled circles) and H*^ 
(open squares), as extracted using Eq. plotted as a function 
of sheet density squared. The best fit of Eq. J3 to H* s (dotted 
curve) is used to extract 7, 01 and 02- Alternatively, the best 
linear fit to H* a (solid line) is used to extract 7. 



gives, within error bars, the same results as fitting over 
the entire measured range of B. Modification of the com- 
mutator [k + 2eA* s , r] by A* s induces spin flipping terms 
~ 7& 3 /4 in the transformed Hamiltonian H*. The cor- 
responding H* D — j^7r 2 m* 2 "/ 2 n 2 /eh, using its expression 
in the diffusive regime. 

We now turn to a discussion of the experiment. 
Samples on three separate heterostructure materials all 
showed qualitatively similar behavior. The sample for 
which data is presented consists of a GaAs/AlGaAs het- 
erostructure grown in the [001] direction with double 8- 
doping layers set back 143 A and 161 A from the 2DEG 
and a total distance of 349 A from the surface to the 
2DEG. A 200 /im wide Hall bar with 700 ^m between 
voltage probes was patterned by wet etching. A litho- 
graphically defined Cr/Au top gate was used to con- 
trol density and mobility in the Hall bar over the range 
n = 1.4-7.0 xlO 15 mr 2 and \i = 3.6-31 m 2 /Vs. Mea- 
surements were made in a 3 He cryostat at temperature 
T = 300 rriK using ac lock-in techniques with bias cur- 
rents ranging from 50 to 500 nA. Figure 1(a) shows 
the longitudinal magnetoconductance as a function of 
V g . A crossover from pure WL (Fig. 1(a), inset) at 
V g = -240 mV to essentially pure AL at V g = +250 mV 
is observed. This crossover demonstrates that a gate can 
be used to control SO over a wide range, as pure WL 
corresponds to negligible SO rotations within the phase 
coherence length L v , while AL corresponds to spin ro- 
tations > 2ir. The solid curves in Fig. 1(a) are fits of 
Eq. J2J) with three free parameters, H$, H* a , and H* s . 
H tr is fixed at each gate voltage by measured values of 
density and mobility. Figure 2 shows extracted param- 
eters H* Q and H* s as a function of n 2 . H* Q is well de- 
scribed by the predicted linear dependence on n 2 , with a 
best fit (Fig. 2, solid line) giving 7 = 31 ± 3 eVA 3 with 
zero y- intercept (see Eq. (tTc|) L The density dependence 
of H* s is well described by Eq. J3J), (Fig. 2, dotted curve), 
giving fit parameters 7 = 28 ± 4 eVA 3 , ot\ = 4 ± 1 meVA 
and «2 = 5± 1 meVA. In this way, the three SO parame- 
ters ai, a.2, and 7 are separately obtained from transport 
measurements by explicitly making use of the density de- 
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FIG. 3: (a) Magnitudes of isotropic linear Dresselhaus (Qui) and 
Rashba (f2j^) terms, and nonisotropic cubic Dresselhaus (£Id3) 
term as functions of gate voltage, V g , density, n, and mobility, 
fi. Insets show theoretical dependence on momentum direction for 
the three terms, indicating that the linear terms are isotropic, while 
the cubic term has a four-fold symmetry and is highly anisotropic. 
Maximum magnitude (when 4> = (j + t)t) is shown for the 
anisotropic (Od3) term, (b) Angular variation of f2, the magni- 
tude of the total SO precession vector at V g = —150 mV (dotted), 
mV (dashed), and 250 mV (solid), corresponding to densities of 
2.3, 5.0, and 7.0 X 10 15 m -2 respectively. 



pendence of H* s and H* a . Extracted values of H v cor- 
respond to dephasing times in the range t v ~ 0.1-1.0 ns 
at 300 mK, which decrease by more than an order of 
magnitude as temperature is increased to 2.5 K. Within 
the error bars, H* Q and H* s do not depend on temper- 
ature over this temperature range. Figure 3(a) displays 
the magnitudes of the three spin-orbit terms as functions 
of V g , n, and //, determined using Eq. and the ex- 
tracted values of a±, ai, and 7. Plotted are values along 
the [110] direction, 4> = tan~ 1 {k y /k x ) = ?, where 0^3 is 
maximum. The total spin precession rate, f2, is plotted as 
a function of the direction, <f>, of the electron momentum 
in Fig. 3(b). While for most directions is an increasing 
function of density, it is seen to decrease with increasing 
density near <fi = ^ and ^ . The linear Dresselhaus and 
Rashba terms (fim and flu) are of comparable magni- 
tude to each other for all densities and in all directions. 
Near (j> = ^ (j an integer), <C £lm,&R and the 

SO is controlled by the linear terms. For <fi near ( 2 - ? '+ 1 ) 7r , 
the cubic term becomes comparable to or even exceeds 
(at high densities) the linear terms. Depending on cj), the 
linear and cubic terms either add {(f) ~ ? , -^p) or subtract 
{<t> ~ X' "T )■ The extracted values for 7 (31 ±3 eVA 3 us- 



4 



ing H* Q , 28 ± 4 eVA 3 using H* s ) are in good agreement 
with the value 27.5 eVA 3 from band structure calcula- 
tions [2iJ |3jj. Values of f2 are ^3 — 8 times smaller 
than previously measured using Shubnikov-deHaas os- 
cillations |26j , with corresponding theory [26l l37| lying 
roughly between the experimental ranges. We note, how- 
ever, that the values are sample dependent. Estimates for 
«i give values for (k^) that correspond to a wave function 
width of ~ 10 nm in the z direction, which is also reason- 
able. The extracted «2 corresponds to a uniform electric 
field E ~ 10 MV/m, using a 2 = cxpeE and a value of 
«o = 5.33 A 2 from a k • p morlel |3fi[ |37| . Previously 
existing models for WL/AL p. Ha l3q i)rovide fits to the 
data that appear qualitatively reasonable, giving values 
for H so that are ~ 5 times higher than those found using 
Eq. (J2J. However, these fits also lead to the unphysi- 
cal result that t so < r. Such unphysical results are not 
surprising given that, for V g > — 50 mV. th e SO length, 
vf /{&>}, is less than £, while theory |allMl29 assumes dif- 
fusive spin evolution £ <C X so , L v . Wo, note that a theory 
for arbitrarily strong SO coupling |l4| may also be used 
to fit this data by including B via L v . This approach 
yields values for fli>3 and Sim consistent with Eq. to 
within a factor of ~ 3, but does not separate Sim and 
S}r terms. We thank I. Alciner, H. Bruus and S. Stu- 
denikin for illuminating discussions and F. Mancoff for 
device fabrication. This work was supported in part by 
DARPA-QuIST, DARPA-SpinS, ARO-MURI, and NSF- 
NSEC. We also acknowledge support from ONR and NSA 
(Y. L.-G.), NDSEG (J. B. M.) and the Harvard Society 
of Fellows (D.G.-G). Work at UCSB was supported by 
QUEST, an NSF Science and Technology Center. 
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